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Abstract 


In the design of thermal power plant components, safety is the prime consideration. 
Generally components are designed as per the guidelines of ASME Boiler and Pressure 
Vessel Codes. In the design, the life of the plant is estimated based on assumed process 
transients. Because of the conservatism in the design actual life of the components is 
more than estimated life. Recently the issue of remnant life of the components has 
attracted considerable attention in the power industry. 

Aging of the components is one of the prime considerations for safe operation of the 
plant and fatigue is the common aging mechanism. Hence to estimate the remnant life of 
the power plant components, on-line monitoring of fatigue degradation is necessary. 

In the present work feasibility study of using Artificial Neural Networks in conventional 
Finite Element Based on-line fatigue monitoring system is done. Feed Forward Error 
Back-propagation neural network was trained and tested for the data obtained by FEM 
code, which converts transient in steam temperature to thermal stress response of the 
component. Effect of over training on the performance of the network is also studied. 



Nomenclature 


a Thermal coefficient of expansion. 

A Area of the element. 

[B] Matrix relating displacements with strains. 

[D] Constitutive matrix, 

s Strain. 

/ Distributed force per unit volume. 

[K] Global stiffness matrix, 

nj Direction cosine of normal to the surface. 

n n Number of stress cycles of stress level <r n , already being seen by the component. 

[N] Matrix of shape functions. 

N n Fatigue life of the component at stress level an. 
t) Learnung rate constanat. 

n Potential energy. 

Pb Body force. 

Ps Surface force, 

a Stress 

t Surface Traction, 

t Transient build up time. 

AT Change in steam temperature before operator takes corrective measures, 
u Displacement. 

U Strain energy, 

v Poison’s ratio. 

Wj ti Weight connecting i th neuron to j* neuron. 

Input vector. 


x 
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Chapter 1 


Introduction 

Computer simulation of the behavior of complex systems and components, whose 
dynamic representation typically requires the solution of many equations and extensive 
use of closure relations, has become very important in modem design. This way of 
proceeding is used particularly in the nuclear industry where safety rules are rather rigid 
and impose to closely examine any possible situation and mode of operation of the 
system. 

The traditional approach to dynamic simulation is based on the realization of a 
simulation model in three steps: understanding the underlying physics, defining the 
mathematical model and developing a fast and reliable code. This is a challenging process 
that needs to be carried out with care. For example, the underlying physics may contain 
significant non-linearities that render the associated mathematical model quite 
complicated. On the other hand, the simplifications and approximations, which are 
necessarily performed, may lead to the introduction of uncertainties. Moreover, in the 
successive phase of implementation of the physical-mathematical model, issues of 
stability and convergence can complicate the discretization of the differential systems. 
Still, the main problem of numerical simulations is that if the equations of the model are 
strongly coupled and the characteristic times of the involved phenomena differ 
substantantially, the traditional methods of solution may involve long, iterative processes 
which result in simulation times significantly larger than the real time of the transient. 

Artificial Neural Networks (ANN) can be used in this context to somewhat 
alleviate this problem. Their parallel and multi-parametric character, and their great speed 
of computing have permitted to assert themselves as a powerful computational tool in 
many fields of research and application, especially where traditional techniques involve 
long elaboration times or when it is particularly difficult to define the underlying 
physical-mathematical model [10]. 

The advantages that result from the ANN as a simulation tool are found not only 
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necessary, in principle, to completely unravel the physics that governs the plant, since 
ANN is capable of constructing its own internal representation of the system input/output 
relations on the basis of the training set. 

Moreover, it is not necessary to build and solve a system of partial differential 
equations but it is sufficient to have some data covering all the possible situations of 
operation, for the ANN training. The generation of training set can be done either by a 
numerical code or from actual plant operation data. What is important is that the training 
data cover all the working space that we are going to simulate. 

The objective of the present work is to study the feasibility of neural network for 
the on-line computation of transient temperature and thermal stresses in a thick pipe. The 
ANN needs training data to train the network which is obtained either by experimentation 
or by numerical simulation. In the present work Finite Element Method (FEM) has been 
used as a numerical technique to obtain the temperature-stress response of the structure, 
which will be used to train and test the network. 

1.1 Problem Formulation 

The Fig. 1.1 shows the schematic of a thick pipe for which temperature and stress 
response of pipe material is predicted using adaptive time delay neural network. Steam at 
constant pressure and high temperature flows through it with constant flow rate. Transient 
Stress and Temperature response of the pipe material is to be simulated for step change in 
temperature. Following steps are to be carried out in the present project: 

Step 1 

Generate training and testing data using FEM, which converts plant transients to 
temperature and stress response of the pipe material. 

Step 2 

Construct and train the network with MATLAB’s NN TOOLBOX for above data and 
optimize network topology. 

Step 3 

Test the above network for new input data not seen during training. 

Step 4 

Predict the remaining life of the component to avoid disaster. 
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Fig. 1.1.1 Schematic of a thick pipe for which stress response of pipe material is predicted 

1.2 Thesis Layout 

• Chapter 2 describes detailed problem description and basics of FEM used to solve 
solid mechanics problems. 

• Chapter 3 gives details neural networks and training and testing data generation 
using FEM. 

• Chapter 4 gives simulation results. 

• Chapter 5 summarizes the conclusion and suggests the scope for future work. 




Chapter 2 


Problem Description 

Main working fluid of thermal power plants is steam, which flows through thick pipes. 
Thermal power plant operates with several plant parameters like temperature, pressure 
and flow rate of the steam and plant is supposed to run at their constant values. But in 
plant operation it’s very difficult to have a control over them. Malfunctioning of any 
plant component results in the fluctuation of these plant parameters. When operator 
senses fluctuation in plant parameter, he takes some corrective action to bring it to their 
normal plant operation values. In the present work the fluctuation in the steam 
temperature is taken into consideration. Fluctuation in the steam temperature will result 
in the cyclic fluctuation in the thermal stress in the pipe, which in turn will result into the 
fatigue of the material. Fatigue is the main aging i.e. degradation of material with time, 
mechanism and hence, common cause for most of the failures in power plant 
components. This effect is the most dominant when the plant approaches the end of the 
designed life. Hence to estimate the remnant life of the power plant components, on-line 
monitoring of fatigue life is necessary [6]. Following are the common methodologies 
used for fatigue monitoring: 

1 . Analytical solution and/or empirical equations as per the standard codes. 

2. Green Function Technique (GFT). 

3. Finite Element Methods. 

On-Line Fatigue Monitoring System based on FEM has been successfully developed and 
commissioned by Reactor Safety Division, BARC, Mumbai. That converts plant 
transients to Temperature-Stress response of the structure using FEM. Computation with 
FEM for On-Line Monitoring of fatigue is computationally expensive and sometimes the 
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time taken by FEM module to convert plant transients to temperature-stress response of 
the structure is more than the duration of the transient itself [7]. To study the feasibility of 
replacing the FEM Module with Artificial Neural Network Module is the main objective 
of the present work. Advantages of such a scheme are as follows: 

1 . ANN can be trained off-line so on-line computational time will be very less in 
comparison with that of FEM. 

2. For better accuracy FEM needs very fine mesh that is analyzed at every node. 
This process creates large matrices and requires large memory space of 
computing machine. Also even a single data corruption can cause trouble to the 
process and ultimately accuracy will suffer. While what ANN needs is a storage 
of weights that is created during training phase. 

2.1 Analysis Procedure 

In the present project On-Line monitoring of a thick pipe as described in Fig 1.1.1 
through which steam flows at high temperature is done at two critical points viz. innei 
radius and outer radius. For on-line monitoring quantities of interest are: maximum stress, 
minimum stress, their difference and the time instant when these occur at both the critical 
points. 

Problem can be divided in three parts: Generation of training and testing data for 
ANN by FEM for various boundary conditions (input to the NN are difference in 
temperature of the steam from normal operating temperature and duration of the transient 
while output of the NN are maximum stress, minimum stress, their difference and the 
time instant when these occur), Training of the Neural Network (code available in 
MATLAB’s NN TOOL BOX) and Testing the Network for the input values not seen in 
training. In the next section FEM methodology is reviewed and in next chapter it’s use in 
training and testing data generation is described. 
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2.2 Finite Element Method 

Finite Element Method (FEM) is a numerical technique to obtain numerical solution to a 
wide variety of engineering problems. Though simple idealized problems have exact 
solutions, in general, the real engineering problems are too complex to permit an exact 
solution. They demand quantitative predictions. They have physical parameters that need 
constant study and are not amenable for exact solution. In such real problems ‘numerical 
approximation is the only resource to obtain quantitative predictions. At present several 
numerical methods are available, namely, finite difference, finite elements, finite 
volumes. Power series, spectral methods, boundary solutions etc. FEM is extensively 
used in various fields such as structural mechanics, solid mechanics, thermodynamics, 
electromagnatics etc. The applications include the design of automobiles, airframes, 
high-rise buildings, spacecrafts, heat engines, electric motors, bearings etc. 

The FEM may be described as a numerical method in which the unknown 
function in the problem domain is approximated by piecewise defined functions. That is, 
the domain is divided into several sub-domains of simple shapes and function is 
approximated separately in each sub-domain by simple polynomial functions. The 
governing equations are derived in a overall sense as a summation of elemental 
equations. The sub-domains are called ‘elements’ and they are interconnected at points 
known as ‘nodes’. 

FEM has its maximum application in structural mechanics and the solid 
mechanics. Solid mechanics deals with the general formulation of deformable bodies 
subjected to forces and displacements, stresses and strains. In following sections basic 
solid mechanics concepts like equilibrium equations, compatibility condition, constitutive 
relations and general FEM formulation for solid mechanics problems are reviewed. 

2.3 Basic Solid Mechanics Concepts 

A three-dimensional body occupying a volume V and bounded by surface S is shown in 
Fig. 2.1. Points in the body are located by x, y, z coordinates. The boundary is 
constrained on some region, where displacements are specified. On part of the boundary, 
distributed force per unit area t , also called traction, is applied. Under the force body 
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deforms. The deformation of a point P (x,y,z) is given by three components of its 
displacements: 

u = [u,v,w] T (2.1) 



Fig. 2.1 General three-dimensional body 


The distributed force per unit volume, for example, the weight per unit volume, is the 
vector / given by: 


f = [xj,z] T 


( 2 . 2 ) 


The surface traction t may be given by its components at points on the surface as: 

t =[/ ,t ,t ] T 

x y z (2.3) 


Stress a acting at a point can be represented by a 3x3 symmetric matrix: 


® [*3x5 ^xy> ^xz? ^yz? ^xz? v yz? 


Vz,] 


(2.4) 


Consideration of force balance in x-, y- and z- directions in bulk V give equilibrium 
equations in x-, y- and z- direction respectively. 


d <rx ( d Jxy [ dr xz 
dx dy dz 


+ X = 0 


(2.5) 
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dx dy 


d _ryz 

dz 


+ Y = 0 


( 2 . 6 ) 


^Ixz+?l21 + ^£l. + Z=0 

dx dy dz (2.7) 


Equilibrium equations on surface S give traction boundary conditions: 


o' x n x + T x y n y + ? xz n z = t x 

(2.8) 

Txy n X + Oyny + Ty Z n z = ty 

(2.9) 

T xz n x + T yz n y + a 'z n z = tz 

(2.10) 


where n x , riy and n z are directional cosines of the external normal on the traction 
boundary r q • 


Strain-Displacement Relations 

Strain s at a point, say P, can be represented by a 3x3 matrix as: 

8 = [Sx> Yxy> Yxzj Yxys £y, Yyz; Yxz; Yyzj S z ] (2.11) 

For small deformations strain s at P can be written in terms of displacement u at P as: 


r du du dv du dw du dy dy dv dw du dw dv dw 
= [ — , — + — , — + — ; — — + — ; — + — + — + — , — ] 
dx dy dx dz dx dy dx dy dz dy dz dx dz dy dz 


( 2 . 12 ) 


Compatibility Condition 

These are the additional conditions based on the fact that the displacements are 
continuous and single valued at every point in the body. From this property compatibility 
condition can be written as: 
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d 2 gy t d^£x _ 
d dy ^ dxdy 


(2.13) 


Constitutive Relations 

For linear elastic materials, the stress-strain relations come from generalized Hooke’s 
law. For isotropic materials, the two material properties are Young’s Modulus of 
Elasticity and Poisson’s ratio. For linear elastic isotropic material stress-strain relation is 
given as: 

"l-i/ v v 0 0 0 I 

v l- V V 0 0 0 £x 

v v \-v o o o £ y 

£ ooo — — o o \ £z 

(1 + v)(l - 2v) 2 >> 

0 0 0 0 0 

oooo o v *y 

L 2 J 




Or, 

M-MM 

where [d] is constitutive matrix. 


2.4 Basic FEM Formulation for Solid Mechanics 
Problems 

There are two commonly used FEM formulations for solid mechanics problems -Galerkin 
Weighted Residual Formulation and Variational approach. In present project variational 
approach has been used hence, will be reviewed here. 

2.4.1 Variational Approach 

In structural mechanics there exist several variational methods such as minimum 
potential energy, maximum complementary energy, stationary principle of Hellinger- 


o 
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Reisner etc. The principle of minimum potential energy is a very popular variational 
principle. That is, out of all displacement fields satisfying displacement boundary 
conditions, the displacement field, that minimizes the potential energy is the one that 
satisfies the governing equilibrium equations. The potential energy II involves the 
internal strain energy of the body and the potential due to external forces. That is, 


n = U-W (2.15) 

where U is the total strain energy of the body and -W is the potential due to external 
forces. Since strain energy density at a point is given by 

dU = —( (J x£x^ m cr y S y^~ <y Z £z J( ~ T yzY yz~^ ^ zxY zx~^ ^xyY X y ) ( 2 - 1 6 ) 

Potential due to external forces is given by 

W = j (Xu + Yv + Zw)dv + | ( t x u + t y v + t.w)dr ^ 

v r q \ ■ J 


Hence, Potential Energy of the body can be given as: 

n = \ (<J X Sx + <Ty £y + CTz Sz + Tyz V yz + T a Y a + Txy YxyJ^V 

A v 

-\(Xu + Yv + Zw)dv- j(t^i + t y v+ J, w)df 

v r 

V 


(2.18) 


Where integrations are carried over the whole volume of the body and on the surface 
where tractions are prescribed. This can expressed concisely using matrix notation as, 

n = -j{x} T {u}dv - J {I } T {u}dr (2 _ j 9) 

v v r q 

In finite element analysis the P.E. of the whole body is expressed as the summation of PE 
over all the elements as 


n= 

<2 = 1 

Potential Energy for one element can be calculated as 

= { \{e} T {o}dv- J {x} T {u}dv- j{t} T {u}dr 


v (e) 


v (e) 


fc 


( 2 . 20 ) 


( 2 . 21 ) 
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or, 


n (e) = L l{e) T Me}- I {x} T {u}dv- J {<} T {«}dr 


,(e) 


/ e ) 


r a 


( 2 . 22 ) 


y\ V' ‘ ' 1 q 

where [D] is constitutive matrix which related stress tensor to strain tensor as given in 
equation (2.10 a) 

Above can be written in terms of displacements as 

sT, 


n< e >= IbY'*' >FM* 

v (‘) v (e> 

- f {uY e ^ [n]^7cIs 


(2.23) 


where [2?] relates displacements with strain and [jV] is the matrix of shape functions. 
Element displacement vector is related with nodal displacement by shape functions as 


(2.24) 


and v (e) is the domain of the element. The expression is known as 

v («> 

T 

elemental stiffness matrix and is represented as , and J is known 

r a 


r -s 

u 


~ N J 

0 

0 ' 

' 

U J 

] v 

► = 2 

0 

N J 

0 

V J 

w 

j 

0 

0 

N J_ 

W j 

- J J 


as elemental force vector p^) . Hence, 


(2.25) 


above elemental matrices can be summed to Global Matrices by the process of assembly 

n =~{uf MM- {uf{Pb\- {uf {p s } 


(2.26) 
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where [AT] is Global stiffness Matrix i.e. matrix obtained by assembling elemental 

stiffness matrices, and 

{«} is global nodal displacement vector. 

As per minimum PE principle, n has to be minimized w.r.t. {«}. 

Hence, considering 

5n fn\ 

4T {0} (2 ' 27) 

gives 

W«}-W}-W}={0} (2.28) 

Or. 

M»}=M (2.29) 

{p} includes concentrated load, if there are any. 

In present problem geometry and load both are axisymmetric, so problem can be 
solved with axisymmetric methods of FEM. Element selected for the analysis is 
triangular ring type element. Next subsection reviews the FEM analysis of axisymmetric 
solids with triangular ring type element. 

2.5 FEM Analysis for Axisymmetric Solids 

There are many practical applications where the bodies are symmetric about one axis of 
symmetry such as, cooling tower, engine cylinder etc. Revolving a plane cross-section 
about the axis of symmetry can generate such bodies. If the forces acting on the solid are 
also axisymmetric then the solution i.e. displacement, strains, stresses etc.) will also be 
axisymmetric. In these cases, geometry, material, loading, boundary conditions are all 
independent of angular coordinate 0 leading to tangential displacement, v=0, and the 
radial and axial displacements u and w are functions of r and z only. Though these 
problems are three-dimensional, it is enough to analyze the radial cross-section (rz-plane) 
only. 

u = u(r,z), v = 0, w = w(r,z) (2.30) 


t 
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Therefore, for axisymmetric problems the strain tensor in polar cylindrical coordinates 
can be written as: 


{,} = 


er Yrd Yrz 
YrQ £0 YQz 
Yrz YQz £z 


du 

dr 

0 

du dw 
dz dr 


0 

u 

r 

0 


du dw 
dz dr 
0 

dw 

dz 


(2.31) 


Its obvious from Hooke’s Law that there will be only four non zero stress components 
and p Td = = 0 . Hence non-zero stress components are related to non-zero strain 


components as: 


(Jr 


'l-v 

V 

V 

0 

Sr 

(Je 

E 

V 

l-v 

V 

0 

Se 

•s 

o-z 

" ~ (l + v)(l-2v) 

V 

V 

1 - 1/ 

0 

1 -2v 

2 . 

Sz 

„ T rz , 


0 

0 

0 

y rz„ 


(2.32) 


In the FE formulation of axisymmetric problems, the element is also axisymmetric ring 
type that is generated by revolving a plane element in rz-plane about z-axis. If the 
axisymmetric element has triangular cross-section in rz-plane and i j and k represent the 
nodes of the triangle. Then displacement components can be written in terms of nodal 


displacements as 

r " 
uj 

Wj 


|M(r,z)l _ 

>/ 

0 

Nj 0 Nk 

0 ' 

[w(/-,z)J 

0 

Ni 

o 

o 

Nk_ 


u ki 

. w k. 


(2.33) 


where Nj , N j and N k are shape functions. And strain components are related to the 
nodal displacements by the relation: 
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b J 

Nj 

0 

bk 

0 " 
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Sr 

se 
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0 

Nk 

0 

Wj 

U J 

Wj 

SO 
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C J 
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. ci 

bi 

c j 

b j 

Cc 
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™k. 


=MM <£) 


(2.34) 


where 


{xjy k -x, \yp 
a/= 2^4 

bl ~ 2 A 

(x k ~xj) 

Ci= 2 A 

ixkyi~xiy k ) 

r C Pk-yi ) 

(xi~xk) 

a r 2 A 

b J~ 2 A 

C J~ 2 A 

(x/y j-xjyj ) 

O' i-yj ) 

(x j - Xj ) 

ak ~ 2 A 

bk ~ 2 A 

ck ~ 2 A 

here (x/.J'/), (x j,y j] 

1 and [xk^k) are coordinates of nodes i, j and k respectively 


and A is the area of the triangular section of the ring type axisymmetric element in rz- 
plane. 

Now elemental stiffness matrix can be calculated as usual by 


[Kp= l[B] T [D\B]dv 

v (0 


(2.35) 


Since quantities are independent of 0, integration on 0 from 0 to 2n can be performed. 
Thus 


[K\( e > = \\[Bf[D\B]2nrdrdz 

Q (e > 


(2.36) 


As the element becomes smaller and smaller, the variation of r and z becomes very small 
compared to the value of r, and hence [5] can be considered to be constant within the 
element. As first approximation stiffness matrix can be determined by considering the 
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integrand to be constant equal to its value at the centroid of the triangle and multiply it 
with area of triangle. 

Hence stiffness matrix is 

[k]^ = [b} T [Df B}lxFA (2.37) 


T 

where [B]‘ is the matrix [5] evaluated at the centroid and A is the area of triangle 
having nodes ij,k. 

In axisymmetric problems forces are also distributed axisymmetrically. Load 
vector due to body force Zj can be calculated using natural coordinates. For 


triangular ring type element elemental load vector for body forces is given by: 



(2 rj + r j + rk)R 
(2 rj + r j + rk)Z 
rtA (ri + 2rj + rk)R 
12 (ri + 2rj + rk)Z 
(ri + r j + 2rk)R 
(r/ + r j + 2rk)Z 


(2.38) 


T 

Load vector due to traction {j r T z } and assuming ij-side of the element coincides 

with traction boundary can be given as: 


f(2 n + r j)Tr 



(2r/ + r j)T z 
(rj + 2r j)Tr > 

( ri + 2r j)T z 


0 

0 


(2.39) 


where Sjj is the length of side ij. 

Using FEM, thermal problems can also be solved following similar steps as that 
of solving solid mechanics problems. In thermal problems having convective boundary 
condition element stiffness matrix gets added with another matrix known as element 
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stiffness matrix due to convective boundary condition for axisymmetric 

triangular ring type element having ij side on the convective boundary, stiffness matrix 
due to convective boundary condition is given by: 


[K}ft = 2xr si] h 


I I 0 

3 6 

- > 0 
6 3 

0 0 0 


( 2 . 40 ) 


where h is convective heat transfer coefficient and r is the radial distance of 
convective boundary. 

Thermal load vector for convective heat transfer problem is given by: 


(2ri + r j)h(t>^ 
0 



0 

0 


here Sij is the length of side ij which is on the convective boundary, tf)^ is the 

temperature of the fluid and h is the convective heat transfer coefficient between fluid 
and surface. 

Next section reviews the application of above theory to calculate thermal stresses in 
axisymmetric solids. 

2.6 FEM Analysis for Thermal Stresses in 
Axisymmetric Solids 

Some times an elastic solid is subjected to deformation due to effects like temperature 
change or moisture absorption etc. For example, the strain in a 1-D rod due to a 
temperature change AT can be given as: 
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£ Q = aAT 


(2.42) 


where a is the thermal expansion coefficient, and AT is the temperature change. For 
triangular axisymmetric element thermal strain is given by: 

ko}={aAr aAT aAT 0} r (2.43) 

If the body is free to expand initial strain will develop due to temperature change and if it 
is constraint to expand initial stresses will develop. A combined expression for state of 
stress considering both initial stress and initial strain can be given as: 

M = MM" feoS+W (2-44) 


where {cro} and {sq} are initial stress and strain respectively. 

Element stiffness matrix [^( e )J can be calculated as usual by: 

M (e, = J[ Bf[D\B\h 

„(«) 

Following the potential energy formulation given in previous section one can get thermal 
load vector as: 


(2.45) 




For very fine mesh thermal load vector can be given by: 
Ni 


(2.46) 




2ri\ 


bi 

0 

b j 

0 

bk 

0 


r 

0 

Nj_ 

r 

0 

Nk 


0 ci 

ci bi 
0 


c j 

J u j 


C i b 


0 Ck 
0 Ck bk\ 


"l — V 

V 

V 

0 


a 

V 

l-v 

V 

0 

E 

a 

V 

V 

l-v 

0 

l-2v 

(1 + v)(l - 2v) 

> 

a 

0 

0 

0 



0 




2 . 




A TrA (2.47) 


at centroid Nj = Nj = Nk = T- 
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In the present project temperature and stress response of the pipe material through 
which steam is flowing at constant pressure with constant flow rate, is simulated for 
linear change in the temperature of the steam. Next section reviews the transient thermal 
analysis that is necessary to solve the problem [1,2]. 

2.7 Transient Thermal Analysis 

The difference between structural and nonstructural dynamic problems is in the order of 
governing differential equation. Compared to structural analysis these problems deal only 
with the first derivative in time. Governing differential equation for heat transfer problem 
is given as [1]: 

MM+[C]W-{F} (2.48) 

Here (j> is unknown temperature, F is force vector and for very fine mesh [c] can be 
given as: 

fcl=s£ 


2 1 1 
1 2 1 
1 1 2 


(2.49) 


where c and p are specific heat and density of the material respectively. A is the area of 
the element and F is the radial distance of the centroid of the element. 

There are many numerical schemes to solve above equations. In the present project 
Neumark method is used which is unconditionally stable. Neumark method can be 
summarized in following steps: 

1. Form the [AT] ^ and [C]^ matrices. 

2. Calculate (f> and 0 from the initial condition i.e. from initial steady state condition 
prevailing before transient. 

3. Select time step At and parameter a and 8 and calculate following integration 
constants: 

1 S 1 

a 0 = —~2 a ^~T~ a2=~— 

aAt z aAt aAt 


18 
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fl3 = - — 1 
j 2a 



A t(S 

<*5 = T — 
2 ya 


\ 

2 

J 


a(j — A/(l — 8) a~i = SAt 

4. Form effective stiffness matrix [£] = [k\ + a\[C]. 

1 . Calculate the effective load at time t + At 

*' + ^ = *' + A ' + [cW'+a 4 «5' ) 

2. Solve for temperature at time t + At: 

^ +A/ }=[i]- V +Ar 

In the present project + ^ at t=0 is the force vector in the steady state i.e. 
before transient starts.[l] 


2.8 Fatigue Monitoring 

It has been recognized that a metal subjected to a repetitive or fluctuating stress will fail 
at a stress much lower than that required causing failure on a single application of load. 
Failures occurring under conditions of dynamic loading are called fatigue failures, 
because these failures occur only after a considerable period of service. This effect is the 
most dominant when the component approaches the end of the designed life. To avoid 
such failures on-line monitoring of critical power plant component is done. In the present 
work Miner Rule is used to calculate the fatigue usage of the component, which relates 
the cumulative damage under cyclic stressing to net work absorbed by the component. 
According to this rule the number of stress cycles applied, expressed as a percent of the 
number to failure at the given stress level, would be the proportion of useful life 
expended. When the total damage defined by this concept reaches 100%, the component 
should fail. 

At any time, say t, if a component has seen m, n 2 , n 3 n n number of stress cycles 

with maximum stress cri, ctz, < 73 , a n respectively and if corresponding to these stress 


IQ 
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levels the fatigue life of the component is Ni, N 2 , N 3 _ N n , then fatigue usage of the 

component at time t is 


nl_ + n2_ + n3_ + m_ 
Nl N2 N3 Nn 


(2.49) 


Due to the statistical nature of the fatigue it is recommended that when fatigue 
usage of the component reaches 0.9, component should be replaced [3]. 

In the present work it is assumed that electronic sensors acquiesce steam 
temperature round the clock and prepares a data file temp.txt that stores all the deviations 
in the steam temperature and duration of the deviation. At 12 O’clock in the night every 
day one code get executed which reads all the entries in temp.txt and calculates the 
fatigue usage of the components as shown in section 4.4. Same code writes fatigue usage 
of the component, date wise in file usage.txt. 


on 


Chapter 3 


Artificial Neural Network 

3.1 Introduction 

The field of Artificial Neural Networks (ANN) was bom in the 1940’s when McCulloch 
and Pitts proposed a computational model based on a simple neuron - like logical 
element. Donald Hebb described a learning rule for adapting the connection strengths of 
these artificial neurons [12]. The Hebb rule or delta rule became the basis of almost all 
ANN research that was to follow. After a period of relative inactivity from early 50’s to 
late 70’s a resurgence of interest came in early 80’s, with the work of Hopfield and his 
associates as well as Rumelhart, Hinton and Williams who developed a generalization of 
the delta rule which made possible the training of powerful multi-layered learning 
networks [12]. 

In the recent years ANNs have been applied to many problems, including pattern 
recognition and classification, associative memory, functional optimization, phoneme 
recognition, word spotting, talker identification, language recognition, signal separation, 
time series prediction and combinatorial optimization. The decade has also shown much 
progress in the development of new learning algorithms, structured neural networks and 
unsupervised learning. The next section reviews the Error Back-propagation Algorithm, 
which is the most commonly used algorithm. 

3.2 Error Back-propagation Algorithm 

Network used for back propagation algorithm is shown in Fig. 3.1. The network functions 
as follows: Each neuron receives a signal from the neurons in the previous layer, and 
each of those signals is multiplied by a separate weight value. The weighted inputs are 
summed, and passed through a limiting function, which scales the output to a fixed range 
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of values. The output of the limiter is then broadcast to all of the neurons in the next 
layer. So, to use the network to solve a problem, the input values are applied to the input 
neurons of the first layer, the signals are allowed to propagate through the network, and 
the output values are calculated. 



Fig. 3.1 Network Topology for Back-propagation algorithm 

Since the real uniqueness or 'intelligence' of the network exists in the values of the 
weights between neurons, we need a method of adjusting the weights to solve a particular 
problem. For this type of network, the most common learning algorithm is called Back 
Propagation (BP). A BP network learns by example, that is, we must provide a learning 
set that consists of some input examples and the known-correct output for each case. So, 
we use these input-output examples to show the network what type of behavior is 
expected, and the BP algorithm allows the network to adapt. 

In BP algorithm essentially, the idea is to perform a gradient descend on the error 
of the network: 


11 
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1 N 0 

E= 2 S (3-1) 

7-1 

where N is the total no. of output nodes, tj is the desired target value and yj denotes net’s 
response to an input vector or pattern. The weights are adjusted to reduce the difference 
of net’s response and the desired target value according to the delta rule by the following 
formula: 


A w j } i = rtf '(a j)(t j-y j) xj ( 3 . 2 ) 

where f\aj) represents how quickly activation changes due to excitation of input 
and 77 represents learning rate. In general, the delta rule may be written as: 

j,i~V 8 j xi (3.3) 

for output nodes: 

8j = f\a j )(tj-y j ) (3.4) 

The errors of output nodes than further propagate backwards to the units in the 
hidden layer(s) and contribute to the change in weights. Consider a hidden node k and let 
/ £ be the set of nodes which node k connects to ( /£ is the fan-out of k), therefore taking 

into account the contribution to the error from all of these nodes in /£ . Hence error 
signal for hidden unit becomes: 

8k~f'(.cik) Z Sjwjk c\ 

jeJk (3 ' 5) 

Above weight updation continues till the error reaches below desired value .[5,12]. 

Simple feed forward networks using BP algorithm can be trained to accomplish 

pattern recognition tasks with complex non-linear boundaries, they are limited to 

processing static patterns - patterns those are fixed rather than temporal in nature. To 

overcome this limitation time delay neural networks and adaptive time delay neural 

networks are proposed [4,8, 11]. 
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3.3 Training and Testing Data Generation 

3.3.1 Training Data Generation 

Training data is generated using FEM code, which converts transient in the steam 
temperature to stress response of the pipe material. FEM code is executed for following 
linear changes in the steam temperature taking over time span of 2 to 10 minutes from 
normal operating temperature i.e. 400 °C: 360 °C, 380 °C, 420 °C, 440 °C, 460 °C. When 
operator comes to know about deviation in the steam temperature he takes some 
corrective actions and steam temperature returns to its normal operating value i.e. 400 °C. 
It is assumed that the temperature of the steam returns linearly and takes same time as 
taken in the deviation from the normal operating temperature. Fig. 3.3. 1.1 shows the 
steam temperature profile during transient. In this figure AT is the change in the steam 
temperature from normal operating temperature before operator takes corrective 
measures to bring the temperature to the normal operating value. In training data 
generation AT = 20°C, 40°C, 60 °C, -20 °C, -40 °C. / is transient buildup time for the 
steam temperature and its values are taken as 2 to 8 minutes for training data generation. 
For each case the value of maximum and minimum principal stresses, and their times (in 
seconds) are noted down at inner and outer radii of the pipe. These values are given in for 
table 3.3. 1.1. 



Fig. 3.3. 1.1 Steam temperature profile during transient 
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Change Time 
(degC) (min) 


Table 3.3.1 .1 Training Data Generated by FEM 


Inner radius Outer radius 


Max Stress I Min stress I Max Stress I Min stress 



- 1.5177 


- 1.4849 


1.4849 

1.4849 

1.4849 

1.4849 


1 

- 1.6132 

262 

- 1.4518 

73 

- 1 . 5324 . 




- 1.4489 


- 1.4450 

199 

BE IMWmr/.yjm 

- 1.4396 

326 

- 1.4375 

390 

- 1.4352 

454 

- 1.4338 

517 

- 1.4323 

581 

- 1.4308 

644 

- 1.4130 

136 

- 1.4052 

199 

- 1.3993 

262 

- 1.3944 

326 


- 1.4849 


■II 

BgfclfJl 

383 

- 

1.4 


537 


1.4535 


1.4535 


1.4534 


1.4533 


1.4533 



- 1.4849 


- 1.4849 


- 1.4849 


- 1.3921 


- 1.4548 


- 1.4543 


- 1.4540 


- 1.4538 
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3.3.2 Sample Plots for Training 

Following are some of the sample plots obtained by FEM showing the transient 
maximum principal thermal stress at inner and outer radius of the pipe: 



time step time step 

Fig,. 3.3. 2.1 Max. Principal Thermal Stress Fig 3.3.2.2 Max. Principal Thermal 

at inner radius for+20°C change in steam Stress at outer radius for +20°C change 

temperature (for 2 min transient) in steam temperature (for 2 min transient). 



time step 


Fig3.3.2.3 Max. Principal Thermal Stress 
at inner radius for+40°C change in steam 
temperature (for 4 min transient) 



time step 


Fig 3.3 .2.4 Max. Principal Thermal 
Stress at outer radius for +40°C change 
in steam temperature (for 4 min transient). 
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lime step 


Fig 3. 3. 2.5 Max. Principal Thermal Stress 
at inner radius for+60°C change in steam 
temperature (for 6 min transient) 



feme step 


Fig 3. 3 .2.6 Max. Principal Thermal 
Stress at outer radius for +60°C change 
in steam temperature (for 6 min transient). 



Urn* step 

Fig3.3.2.7 Max. Principal Thermal Stress 
at inner radius for -20°C change in steam 
temperature (for 8 min transient) 



time step 


Fig 3. 3 .2. 8 Max. Principal Thermal 
Stress at outer radius for -20°C change 
in steam temperature (for 8 min transient). 
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Fig 3. 3. 2.9 Max. Principal Thermal Stress 
at inner radius for -40°C change in steam 
temperature (for 5 min transient) 

3.3.3 Testing Data Generation 



Fig 3.3.2.10 Max. Principal Thermal 
Stress at outer radius for -40°C change 
in steam temperature (for 5 min transient). 


The same FEM Code that generates training data for the input parameters during training 
also generates testing data. Testing data is generated for AT = 30°C, 50°C, 70 °C, -30 °C, 
-50 °C and t = 3min, 4min, 5min and 6min. These parameters lie out of the parameters 
used for generating training data. Testing data is given in table below in the same way as 
the training data and will be used to compare the simulation results of neural network in 
the next chapter. 


Table 3.3.3. 1 Testing Data Generated by FEM 


Temp 

Change 

(degC) 

Tran 

Time 

(min) 

Inner radius 

Outer radius 

Max Stress 

Min stress 

Max Stress 


Value 

Time 

Value 

Time 

Value 

Time 

Value 

Time 

flPMiai 

(s) 


(8) 


(8) 

lalfliml 

(8) 


3 

-1.4849 

1 


199 

-1.4522 

73 

-1.4854 

383 


4 


1 




73 

-1.4928 

482 

5 


1 

-1.5527 

326 

-1.4526 

73 

-1.4989 

567 


6 

-1.4849 

1 

-1.5590 

390 

-1.4527 

73 

-1.5041 

647 


3 

-1.4849 

1 


199 

-1.4517 

73 

-1.5069 

383 

ten 

4 

-1.4849 

1 

-1.5918 

262 . 


73 


482 

«OU 

IT” 


1 

-1.5980 

326 

-1.4522 

73 

-1.5293 

567 


6 

-1.4849 

1 

-1.6033 

mm 


73 

-1.5381 

647 



































































































Chapter 4 


Results and Discussion 

4.1 Parameters 

ANN with back propagation learning paradigm is used to train the network to match input- 
output relationships using the training cases discussed in the previous chapter. Following 
parameter values are used: 

4.1.1 Initial Weights 

The weights of the network to be trained are typically initialized at small random values. The 
initialization strongly affects the ultimate solution. If all the weights start out with equal 
values, and if the solution requires that the unequal weights be developed, the network may 
not train properly. Unless the network is disturbed by random factors or random character of 
input patterns during training. The choice of initial weights is, however one of the several 
factors affecting the training of the network towards an acceptable error minimum. 

4.1.2 Learning Constant 

The effectiveness and convergence of the error back-propagation learning algorithm depends 
significantly on the values of the weight learning rate constants i.e. 77, . In general, however, 

the optimum value of learning rates depends on the problem being solved and there is no 
single learning constant value suitable for different training cases. 

4.1.3 Momentum Factor 

The purpose of momentum factor is to accelerate the convergence of the training. The factor 
involves supplementing the current weight adjustment with a fraction of the most recent 
weight adjustments. 

4.2 Training Exit Criterion 

Training of the neural network stops if any one of the following conditions is met during 
training: 
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• Error level reaches a pre-specified error value. 

• All specified training epochs get exhausted. 

4.3 Results 

A feed forward back propagation neural network with the following network topology was 
trained with the training data generated by FEM Module and given in the table 3.3. 1.1. 

A three-layered fully connected NN (one input layer, one hidden layer and one output 
layer) was configured with two input neurons taking temperature change and duration of the 
transient as input to the NN, eighteen hidden neurons and eight output neurons giving 
maximum and minimum stress and time (in sec) of their occurrence at both inner and outer 
radii. All the neurons use sigmoidal squashing function. Other network parameters are given 
below: 

Learning Rate for weights = 0.4 
Momentum Factor for weight updation = 0.9 
Maximum no. of epochs = 50,000 
Network is trained in batch mode. 

Fig 4.3.1 to Fig 4.3.3 show variation of mean squared error (MSE) vs. no. of epochs starting 
from different random initial weights, as training proceeds. Its obvious from the error plots 
that even if the starting error is different for all the three cases but in 50,000 epochs training 
network converges to approximately same error level. 

Table 4.3.1 shows the simulation results for the testing data shown in obtained in previous 
chapter. Negative errors show under-prediction. 

In the present work the data is normalized linearly between maximum and minimum 
values of parameters. While predicting time instants for maximum and minimum stresses, 
some of the time values came out to be negative indicating that stresses become maximum or 
minimum before the start of the transient. To avoid such situations following mapping rules 
are applied in addition to normalization of the data: 

• In the prediction of time of occurrence for maximum or minimum stresses at inner 
radius, if predicted time is less than 100 seconds it is mapped at Is. 

• In the prediction of time of occurrence for maximum or minimum stresses at outer 
radius, if predicted time is less than 200 seconds it is mapped at 73s. 


MSE 
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Fig 4.3.1 MSE vs. No. of Epochs for first set of initial weights 


MSE 
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Fig 4.3.2 MSE vs. No. of Epochs for second set of initial weights 
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Temp Tran 
Change Time 
(degC) (min) 


Table 4.3.1 Simulation results for Testing Data 

Inner radius 

Max Stress \ Min stress ] 



Predicted 


Exact 


Value 

(10 8 Pa) 

-1.4832 


-1.4849 


Time 

1 


1 


Exact 


1.4849 



Value 

(10 8 Pa) 

-1.5360 


-1.5446 




















.5937 


1.5844 


0.58 


Time 

205 


199 




Outer 

Max Stress 
Value Time 
(10 8 Pa) (s) 

-1.4535 73 


-1.4522 73 


□ 

82 

EESIill 

199 

BEaai 


Exact 


%Error 


Predicted 


Exact 


%Error 


Predicted 


%Error 


Predicted 


Exac 


%E 


Predicted 


Exact 


%Error 


•1.4809 


•1.4849 


- 0.02 


1.4893 


•1.4849 


-0.29 


1.4874 


1.4849 


0.16 



-0.77 -0.38 


-1.5871 333 


-1.5980 326 


1.09 2.14 


-0.15 


1.4549 


1.4522 


0.18 


1.4538 



-1.4937 


-1.4849 


.59 


1.6135 


1.6346 


.29 


■1.4571 


-1.4516 


.37 



Exact 

-1.4849 

1 

%Error 

-0.48 

0 



Exact 


%Error 



-1.4849 I 1 


- 0.11 0 


-1.4212 211 


-1.4251 199 


0.26 




Predicted 


Exact 


%E 


Predicted 


■1.6507 


1.41 


■1.4853 


■1.4849 


.026 


■1;4847 


•1.4849 


-0.013 


•1.4814 



1.4516 


- 0.02 


1.4442 


1.4521 


-0.54 


1.4170 


1.4208 


-0.26 


radius 

Min stress 
Value Time 
(10 8 Pa) (s) 



0.042 0.622 


-1.4854 

382 

WEBMl 

3.14 

memi 

483 

-1.4928 1 

481 

■mi 

0.41 

-1.5065 

536 

BHH 

566 

-0.33 

-5.3 

-1.5132 

615 

-1.5041 

646 

0.60 

4.78 

-1.5095 

EEf 1 

-1.5069 

EH 

-0.11 

-4.7 

-1.5269 

487 

-1.5192 

482 

0.506 

1.04 

-1.5238 

575 

-1.5293 

567 

-0.35 

1.41 

-1.5325 

648 

-1.5381 

647 

-0.36 

0.15 

-1.5044 

353 

-1.5284 

383 

-1.57 

-7.8 

-1.5248 

494 


482 

EEII 

2.5 

-1.5610 

617 

-1.5598 

567 

WEESM II 

8.81 

-1.5902 

705 

-1.5722 

647 

1.14 

8.96 

BEE3II 

MEM\ 

BE2U9I 

■| 

0.089 

0 

-1.4539 

73 

-1.4538 

73 

ittfll 

m 

-1.4500 

73 
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Exact 

-1.4170 

326 


1 


567 

-1.4536 




%Error 

1.00 

8.89 

-0.23 

0 

0.56 

6.34 

-0.24 

0 



Predicted 

-1.4245 

441 

-1.4812 

1 


708 

-1.4497 

73 


6 

Exact 

-1.4138 

391 

-1.4849 

1 

-1.4021 

647 

-1.4535 

73 



%Error 

0.75 

12.78 

-0.25 

0 

0.49 

mm 

-0.26 

0 



Predicted 

-1.3810 

163 

-1.4855 

1 

-1.3940 

319 

-1.4560 

73 


3 

Exact 

-1.3853 

199 

-1.4849 

1 

-1.3993 

383 

-1.4545 

73 



%Error 

-0.31 

-18.1 

0.04 

0 

-0.37 

-16.1 

0.10 

0 



Predicted 

-1.3678 

279 

-1.4849 

1 

-1.3785 

483 

-1.4557 

73 


4 

Exact 

-1.3779 

262 

-1.4849 

1 

-1.3870 

482 

-1.4542 

73 

-50 


%Error 

-0.73 

6.5 

0 

0 

-0.61 

0.20 

0.10 

0 


Predicted 

-1.3545 

304 

-1.4905 

1 

-1.3679 

521 

-1.4608 

73 


5 

Exact 

-1.3717 

326 

-1.4849 

1 

-1.3769 

567 

-1.4540 

73 



%Error 

-1.15 

6.74 

0.37 

0 

-0.65 

-8.11 

0.46 

0 



Predicted 

-1.3508 

342 

-1.4928 

1 

-1.3645 

585 

-1.4619 

73 


6 

Exact 

-1.3664 

390 

-1.4849 

1 

-1.3681 

647 

-1.4538 

73 



%Error 

0.83 

-14.0 

0.53 

0 

-0.26 

-9.5 

0.55 

0 


4.4 Pattern Memorization by Neural Network 

In neural networks it seems that lower the error level reached in training, the testing results 
will be better. But some times if the error in training reaches below a certain level, network 
looses generalization capability and is said to be over trained. In that case simulation results 
for the pattern used for the training are very accurate but those for the patterns not seen 
during the training are very poor. This is what happened in the present problem when 
network was trained with Levenberg-Marquardt Back-propagation Algorithm. This algorithm 
uses adaptive learning rate and converses very fast and in 2500 epochs only the Mean Square 
Error (MSE) reached to the order of 10' 7 , but with this error level network got over trained. 
Simulation results for the training patterns were very close to the exact values obtained by 
the FEM module but those for the pattern not seen in the training were very poor. Another 
network with same parameters other than the exit error criteria (training gets stopped when 
Mean Square Error (MSE) reaches 0.001) was trained with Levenberg-Marquardt Back- 
propagation Algorithm. Network reached the desired error level in less than 150 epoch 
training. Fig 4.4.1 and fig 4.4.2 gives plots for Mean Square Error (MSE) as training 
progresses for over training and adequate training cases respectively. The simulation results 
for the training patterns and patterns not seen during training with this network are with in 
the acceptable error range. 
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Fig 4.4.1 MSE vs. No. of Epochs for over training of the network 
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4.4.1 Simulation Results for Over Training Case with Training Patterns as 
Input 

Fig 4.4. 1.1 to Fig.4.4.1.4 show some of the simulation results for the patterns seen during 
training. Single hashed bar shows the exact value obtained by FEM while double hashed bar 
shows value predicted by ANN respectively. It can be observed from the figures that the 
network follows the original values nicely. 

Maximum error level in the prediction of various parameters is given in the following table 
4.4. 1.1. Only the magnitude of the error is given not sign. 

Table 4.4. 1.1 Max. error in the prediction of various parameters when training patterns are used as input 

parameters for over training case 


Parameter 

Max. Percentage Error in simulation 

Max. Principal Stress at inner 

radius. 

0.0207 

Time of occurrence of Max. 

Principal Stress at inner radius. 

0.7353 

Min. Principal Stress at inner 

radius. 

0.0377 

Time of occurrence of Min. 

Principal Stress at inner radius. 

0.7353 

Max. Principal Stress at outer 

radius. 

0.0141 

Time of occurrence of Max. 
Principal Stress at outer radius. 

0.6224 

Min. Principal Stress at outer 

radius. 

0.0696 

Time of occurrence of Min. 
Principal Stress at outer radius. 

0.6833 


Principal stresses shown in the following figures are compressive in nature. Here only 
the magnitude of the max. principal stress is given not sign. 
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Fig 4.4. LI Simulation results for maximum of Maximum Principal Stress at inner radius 
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Fig 4.4. 1 .2 Simulation results for time of occurrence of maximum of Maximum Principal Stress at inner radius. 
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Fig 4.4. 1.3 Simulation results for minimum of Maximum Principal Stress at inner radius 


□ Exact Value 
i ■ Predicted Value 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 

Training Case Number 

Fig 4,4 J A Simulation results for time of occurrence of minimum of Maximum Principal Stress at inner 

s 




Results and Discussion 


4.4.2 Simulation Results for Over Training Case with Test Patterns as 
Input 

Fig 4.4.2. 1 to Fig.4.4. 1.4 show some of the simulation results for the test patterns i.e patterns 
not seen during training. Single hashed bar shows the exact value obtained by FEM while 
double hashed bar shows value predicted by ANN respectively. It can be observed from the 
figures that the network’s performance is very poor in this case. 

Maximum error level in the prediction of various parameters is given in the following table 
4.4. 1.1. Only the magnitude of the error is given not sign. 


Table 4.4.2. 1 Max. error in the prediction of various parameters when testing patterns are used as input 

parameters for over training case. 


Parameter 

Max. Percentage Error in simulation 

Max. Principal Stress at inner 

radius. 

20.9 

Time of occurrence of Max. 

Principal Stress at inner radius. 

9.9e3 

Min. Principal Stress at inner 

radius. 

17.41 

Time of occurrence of Min. 

Principal Stress at inner radius. 

7.62e3 

Max. Principal Stress at outer 

radius. 

16.05 

Time of occurrence of Max. 

Principal Stress at outer radius. 

1.92e3 

Min. Principal Stress at outer 

radius. 

24.37 

Time of occurrence of Min. 

Principal Stress at outer radius. 

1.99e3 


Principal stresses shown in the following figures are compressive in nature. Here only 
the magnitude of the max. principal stress is given not sign. 
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Principal stresses shown in the following figures are compressive in nature. Here only 
the magnitude of the max. principal stress is given not sign. 
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Fig 4.4.2. 1 Simulation results for maximum of Maximum Principal Stress at inner radius 
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Fig 4.4 .2.2 Simulation results for time of occurrence of maximum of Maximum Principal Stress at inner radius. 
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Fig 4.4,2,4 Simulation results for time of occurrence of minimum of Maximum Principal Stress at inner radius 

Its clear from the previous sections that when MSE level is of the order of 10' 7 
simulation results for training patterns are very accurate and percentage error is always 
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less than 1 /o, but those for test data are very poor. Hence, at this MSE level network lost its 
generalization capability. Next section gives the simulation results for adequate training case 
i.e. when MSE is of the order of 10' 2 . 

4.3.3 Simulation Results for adequate Training Case with Test Patterns as 
Input 

Fig 4.4.3. 1 to Fig.4.4.3.8 show the simulation results for the test patterns i.e patterns not seen 
during training for adequate training case. Single hashed bar shows the exact value obtained 
by FEM while double hashed bar shows value predicted by ANN respectively. It can be 
observed from the figures that the network’s performance is very poor in this case. 

Maximum error level in the prediction of various parameters is given in the following table 
4.4.3. 1. Only the magnitude of the error is given not sign. 

Table 4.4.3. 1 Max. error in the prediction of various parameters when test patterns are used as input parameters 

for over training case. 


Parameter 

Max. Percentage Error in simulation 

Max. Principal Stress at inner 

radius. 

1.08 

Time of occurrence of Max. 

Principal Stress at inner radius. 

14.57 

Min. Principal Stress at inner 

radius. 

1.42 

Time of occurrence of Min. 

Principal Stress at inner radius. 

13.06 

Max. Principal Stress at outer 

radius. 

1.18 

Time of occurrence of Max. 

Principal Stress at outer radius. 

11.48 

Min. Principal Stress at outer 

radius. 

3.04 

Time of occurrence of Min. 

Principal Stress at outer radius. 

11 
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Principal stresses shown in the following figures are compressive in nature. Here only 
the magnitude of the max. principal stress is given not sign. 


[h Exact Value 
; 1 Pr edicted Value 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

Test Case Number 
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Principal stresses shown in the following figures are compressive in nature. Here only 
the magnitude of the max. principal stress is given not sign. 
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4.4 Fatigue Usage of the Component 

In all the above-simulated cases stress level is below the endurance limit of the material, 
hence these plant transients will have no effect on the fatigue life of the component and 
component can survive for infinite such stress cycles. But for the demonstration purpose let 
us assume that fatigue life corresponding to 20°C, 40°C and 60°C is Nj- 10 6 cycles, jV) = 10 5 
cycles and N 3 = 10 4 cycles respectively. At any particular day if the entries in the data 
acquisition file are: 

Temp Change (AT 0 C) 20 40 60 

No. of Cycles ( n; ) 4,500 4,370 2,860 

It shows that by that day component has seen 104500 transients in the steam temperature of 
intensity 20°C i.e. steam temperature fluctuated from 400 to 420°C, 4370 transients in the 
steam temperature of intensity 40°C and 2860 transients in the steam temperature of intensity 
60°C 

Using equation (2.49) and the above N t and «, values, code calculates the fatigue usage of 
the component to be 


2-^ = 0.3342 

mN, 

Assuming that maximum usage factor of the component before replacement is 0.9, at the end 
of the day remaining fatigue usage is only 0.5608. 



Chapter 5 


Conclusion and Future Scope 

This work was carried out to study the feasibility of using Artificial Neural Network in 
FEM based On-Line Fatigue Monitoring system developed by BARC. Following 
conclusions are drawn from the analysis of various results obtained in previous chapter: 

• Error Back-Propagation Neural Network can be used to predict maximum 
principal stress and its time of occurrence in thick pipe, through which steam is 
flowing at high temperature, for transient change in steam temperature. Maximum 
percentage error in the prediction of various parameters as: 



Max.% error 

Corresponding to Change in 

Steam Temperature (AT°C) 

Max. Princi. Stress at inner radius 

1.15 

-50 

Time of occurrence of Max 

Princi.. Stress at inner radius 

18.09 

-50 

Min. Princi. Stress at inner radius 

1.62 

70 

Time of occurrence of Min 

Princi.. Stress at inner radius 

12.88 

70 

Max. Princi. Stress at outer radius 

0.65 

-50 

Time of occurrence of Max 

Princi.. Stress at outer radius 

16.71 

-50 

Min. Princi. Stress at outer radius 

1.5703 

70 

Time of occurrence of Min 

Princi. .Stress at outer radius 

8.9645 

70 


• Maximum percentage error magnitude is either corresponding to AT (Change in 
steam from normal operating temperature before operator takes a corrective 
measure ) equal to 70°C or -50°C, which lie out of range of AT used for training 
the network. This shows that network is poor in extrapolation. Hence, full possible 
range of input parameters should be used to train the network. 

















• Error is more in prediction of time. Maximum error found is about 1 8%, However 
in the present problem, we are more concerned with the stress level values then 
the time when the maximum stress occur. In predicting maximum stress values, 
the error is less than 2%, which may be well acceptable, considering other sources 
of errors which may arise from measurements etc. In predicting fatigue life 
another important aspect is the number of times maximum stress exceeds the 
fatigue strength according to S-N diagram. Since the developed ANN predicts the 
maximum stress close to that obtained by FEM, reliable prediction of fatigue life 
may be possible with developed ANN. 

• For the present problem if error level during training is reduced beyond a certain 
level, network gets over trained i.e. it starts memorizing training patterns and 
looses generalization capability for input / output relationship. In present problem 
value of MSE of the order of 10‘ 3 is adequate for training LM Back-propagation 
Neural Network. 

5.1 Scope for the Future Work 

Problem solved is essentially temporal in nature. In the present work error back- 
propagation algorithm is used, which gives excellent results in processing static patterns. 
In future temporal pattern processing networks like time -delay neural network or 
adaptive time-delay neural network can be used to solve the problem. 
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Simulation Result Tables 


Table 1, Table 2 and Table 3 give the simulation results for over training case with 
training patterns as input, with test patterns and those of with the test patterns 
respectively. Negative errors show under-prediction of the data. 


Table A.1 Simulation results for training patterns ( over training case ) 


Temp 

Change 

(degC) 

Tran 

Time 

(min) 


Inner radius 

Outer radius 

Max Stress 

Min stress 

Max Stress 

Min stress 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 

+20' 

2 

Predicted 

-148.47 

1 

-152.09 

135 

-145.23 

73 

-146.89 

288 

Exact 

-148.49 

1 

-152.08 

136 

-145.22 

73 

-146.85 

286 

% Error 


0 

0.0066 



0 

0.0272 

0.699 

3 

Predicted 

-148.48 

1 




73 

-147.45 

386 



1 



-145.25 



383 

% Error 

-0.0067 

0 

0.0066 

0 


0 

-0.0068 

0.783 

4 

Predicted 

-148.49 

1 



BElldrl 

73 

-147.93 

479 

Exact 

-148.49 

1 

-152.76 

262 

-145.27 

73 

-147.95 

482 

% Error 

0 

0 

0.0131 

0.38 



-0.0135 

-0.622 

5 

Predicted 

-148.49 

1 





-148.35 

565 

Exact 

-148.49 

1 

-153.01 

326 

-145.28 

73 

-148.36 

568 

% Error 

0 

0 

0.0131 


0 

0 

-0.0067 

-0.528 

6 

Predicted 

-148.48 

1 

-153.23 



73 

-148.71 

647 

Exact 

-148.49 

1 

-153.28 

390 

-145.28 

73 

-148.71 

647 

% Error 

-0.0067 

0 

-0.0326 




0 

0 

7 

Predicted 

-148.49 

1 




73 

-149.05 

723 

Exact 

-148.49 

1 




73 

-149.02 

726 

% Error 



0.0065 

-0.22 


0 

0.0201 

-0.413 

8 

Predicted 

-148.49 


-153.58 

518 

-145.28 

73 

-149.31 

799 



1 

-153.59 

517 

-145.29 

73 

-149.29 

799 

% Error 

0 

0 

-0.0065 

0.19 



0.0134 

0 

9 

Predicted 

-148.49 

1 



-145.29 

73 

-149.54 

872 

Exact 

-148.49 

1 



-145.29 

73 

-149.53 

871 

% Error 


0 

-0.0065 

0.17 

0 

0 

0.0067 

0.115 

10 

Predicted 

-148.50 

1 



-145.31 

73 

-149.73 

942 

Exact 

-148.49 

1 


jilfjSi 



^149.75 

msm 
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Esmmmi 


Exact -148.49 


% Error 



Exact i -148.49 


% Error 0 


-155.69 


-0.0064 


-156.46 


-156.45 


0.0064 



-145.14 


-145.14 


0 


-145.20 


-145.20 


0 




Exact -148.49 



Predicted 

-148.50 

Exact 

-148.49 



-157.53 


145.25 


- 0.0 


-148.45 


-148.40 


0.0337 


-149.61 


-0.0067 


-150.56 


-150.60 


-0.0266 


-151.39 



Predicted -148.50 


Exact -148.49 


% Error 0.0067 


Exac 


% Error 


-157.53 

326 

-145.24 

73 

0 

0 

0.0069 

0 

-157.95 

389 

-145.26 

73 


-0.0063 


-145.25 


0.0069 


-152.11 


0 




Predicted -148.49 


-158.69 


Exact 
% Error 


-148.49 

0 



Predicted 

-148.49 

9 

E 

S39IEESB3 


% Error 

0 



-145.27 

73 

-0.0069 

0 

-145.27 

73 




Predicted 

-148.49 

■ 

Exact 

-148.49 



0.0376 

0 

-159.20 

135 


-145.27 


644 -145.27 


-145.28 


69 


0.0196 


-153.75 


-153.75 


0 


-154.13 


-154.18 


% Error 


Predicted 


Exact 


% Error 


Predicted 


Exac 


% Error 


Predicted 


Exa 


% Error 


Predicted 


Exact 


% Error 


Predicted 

Exact 


0 


-148.50 


-148.49 


0.0067 


-148.50 


0.0067 


-148.50 


-0.0377 


-145.05 


-0.74 -0.0069 


-145.14 






0.0067 


-148.49 


-148.49 


-148.48 

-148.49 


-162.08 326 


-162.06 326 


0.0123 0 


-145.21 


-145.21 



0.0123 


-163.26 454 


45.22 


0 


145.23 



-153.22 


-153.24 


-0.0130 


-154.48 


-154.46 


0.0129 


0 


285 


286 


-0.35 


386 


383 


0.783 


480 


482 


-0.415 


567 


568 


76 


8 


647 


0.155 


-158.33 

453 

-145.25 

73 

-152.74 

725 




55.55 

6 

55.5 

n 

6 


0.0257 


-156.45 




567 


568 


-0.176 


0.154 



52 
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0 | - 0.0122 | 0 




9 Exact I -148.49 1 1 1 - 164.26 | 581 

% Error -0.0067 0 



Exact 



% Error 

-0.0135 

0 

Predicted 

-144.86 

136 

Exact 

-144.89 

136 


148.50 



Exact I -144.50 


% Error 


Predicted -144.22 262 


Exac 


0.0069 


Predicted -143.96 325 








% Error 


Predicted 


Exact 


Predicted 


Exact 


% Error 


Predicted 
Exact 
% Error 


Predicted 


% Error 


Predicted 


Exact 


% Error 


Predicted 


Exact 


0 


-143.75 


-143.75 


-143.52 


-143.52 


0 


-143.38 

-143.38 

0 


-143.23 


-143.08 


-143.08 


-141.32 


-141.30 


-0.31 


390 


390 


0 


454 


454 


0 


518 

517 

0.193 


582 




-148.49 


-00067 


-148.49 


-148.49 


0 


-148.49 


-148.49 


0 


-148.48 



Exact I -139.93 


- 0.0 


-148.49 


-148.49 


0 


-148.50 


-148.49 


-148.49 


-148.49 


0 


-148.49 


-148.49 


0 


-148.49 


-0.0069 

0 

-145.25 

73 

-145.24 

73 

0.0069 

0 

-145.25 

73 

-145.25 

73 

0 

0 

-145.26 

73 


73 

0 

0 



-0.0139 

0 

-143.18 



383 

bwiinitosll 

-142.68 

479 


0.0070 

-0.62 

EEES£IEs9i 

msm 

568 

-0.0140 

-0.18 

-141.34 

647 

-141.34 

647 

0 

0 

-141.46 

725 


0.0071 

-0.14 

-141.33 

799 

-141.33 

799 

0 

0 

BEH.feMI.H1 

-141.09 | 

871 




0.0064 


-157.25 


-157.25 


0 


-157.86 



-0.0696 


58.35 


58.26 


0.0568 


-145.41 


-145.40 


0.0 


45.38 


-145.37 


0.0069 


-145.36 


5.35 


0.0069 


-145.35 


-0.0069 


-145.34 


-145.34 


0 


-145.35 


-145.33 


0.0138 


-145.32 


-145.33 


-0.0069 


-145.32 


-145.33 


-145.33 


5.33 


0 


-145.50 


-145.48 


0.0137 


-145.41 


-145.43 


38 


-145.40 


0 


798 


799 


-0.125 


871 


0 
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942 


-0.106 



53 
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I^E!3&9 


mm 


Kill 


-0.41 


0 



Predicted 

-139.44 

326 


■ 

-139.21 

567 

-145.37 

73 


5 

Exact 

-139.44 

326 

-148.49 

1 

-139.21 

568 

-145.38 

73 



% Error 

0 


0 

0 

0 

-0.17 





Predicted 

-139.01 

389 

-148.49 

1 



-145.37 

73 



Exact 

-139.01 

390 

-148.49 

i 



-145.37 

73 







0 

0 

-0.15 

0 

0 






-148.49 

1 


mm 

-145.38 

73 



Exact 


454 

-148.49 

1 

-137.89 

mm 


73 



% Error 

-0.0072 

0.22 

0 

0 

0 






Predicted 

-138.28 

517 


1 



-145.35 

73 




-138.28 



1 



-145.35 

73 






0 

0 

0 

0 

0 

0 



Predicted 

-137.97 

583 

-148.48 

1 

-136.88 



73 


9 

Exact 

-137.97 

581 

-148.49 

1 



-145.35 

73 



% Error 

0 

0.344 

-0.0067 




0.0069 

0 







1 



-145.34 

73 




-137.68 

644 

-148.49 

1 



-145.35 

73 



% Error 

0 

-0.31 

0.0067 

0 

0 

-0.21 

-0.0069 

0 


Table A.2 Simulation results for the test i.e patterns not seen during training (over training case ) 


Temp 

Change 

(degC) 

Tran 

Time 

(min) 


Inner radius 

Outer radius 

Max Stress 

Min stress 

Max Stress 

Min stress 

Value 

(10 8 Pa) 

Time 

(s) 

Value 

(10 8 Pa) 

m 

Value 

(10 8 Pa) 

Time 

(s) 

Value 

(10 8 Pa) 

Time 

(s) 

+30 


Predicted 

-142.22 

1 

-155.34 

All 

-138.64 

1011 

-143.79 

752 

Exact 

-148.49 

1 

-154.46 

199 

-145.22 

73 

-148.54 

■EE»»1 

%Error 

-4.2225 

0 

0.57 

1 14.573 

-4.531 

1284.93 

-3.198 

96.859 


Predicted 

-143.48 

1 

-155.26 

437 

-139.74 

900 

-144.63 

1002 

Exact 

-148.49 

1 

-154.90 

262 

-145.25 

73 

t 149.28 



-3.374 


0.232 

66.794 

-3.793 

1132.877 

-3.115 

■blillFJ 

5 


-144.50 

1 

-155.20 

- 457 

-140.72 

801 

-145.46 

1217 

Exact 

-148.49 

1 

-155.27 

326 

-145.26 

73 

-149.89 

566 

%Error 

-2.687 

0 

-0.045 

40.184 

-3.125 

997.260 



6 


-145.31 

1 


488 

-141.54 

719 

-146.22 

1412 



1 

-155.90 

391 

-145.27 

73 

-150.41 

B*£Bi 






-2.568 

884.932 

-2.785 

BEBEl 


3 



100 

-158.94 

90 

-147.38 

186 

-154.03 

1030 

Exact 

-148.49 

1 

-158.44 



73 

-150.69 

383 

%Error 

1.111 

9900.0 

0.316 

-54.774 

1.522 

154.794 

2.216 

168.929 

■ 

Predicted 

-149.80 

1 

-159.96 

178 

-147.00 

226 

-155.23 

1223 

Exact 

-148.49 

1 


262 

-145.20 

73 

-151.92 

482 
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0 

0.490 

-32.061 

1.239 




5 

Predicted 

-149.54 

1 


261 


257 

-156.24 


Exact 

-148.49 

1 

-159.80 

326 

-145.22 

73 

-152.93 

■Safi 

%Error 

0.7071 


0.62 

-19.939 

1.019 

252.055 

2.164 


6 



1 




281 



Exact 

-148.49 

1 

-160.33 

391 

-145.24 

73 

-153.81 


%Error 

0.5657 

0 

0.717 

-13.043 

0.839 

284.932 

2.126 

i mm 

+70 

3 



1 

-190.71 

845 

-133.68 

73 

-148.83 

Hrfl 

Exact 

-148.49 

1 

-162.42 

199 

-145.11 

73 


mmm 

%Error 

-11.415 

0 







4 

Predicted 


1 

-187.32 

831 

-133.99 

1476 

-150.44 

wism 


-148.49 

1 




73 


482 

%Error 

-10.108 

0 

14.597 

217.176 

-7.695 

1921.918 



5 

Predicted. 

-135.60 

1 

-183.68 

811 

-134.55 




Exact 

-148.49 

1 

-164.32 

326 

-145.16 

73 

-155.98 


%Error 

-8.681 

0 

11.782 

148.773 

- 7.309 

1836.986 



6 

Predicted 

-137.70 

1 

-180.23 

788 

-135.40 

1324 

-153.17 

sis 

Exact 

-148.49 

1 

-165.07 

391 

-145.21 

73 

-157.22 


%Error 

- 7.266 

0 

9.184 

101.535 

-6.756 

1713.699 



-30 

3 

Predicted 

-171.74 

1 

-144.92 

760 

-162.87 

73 

-180.61 


Exact 

-142.51 

199 

-148.49 

1 


383 



%Error 

20.51 

-99.50 

-2.404 

75900 

14.633 

-80.94 


WBBEk 

4 

Predicted 

-171.77 

1 

-144.88 

763 

-162.87 

73 

-180.81 


Exact 


262 

-148.49 

1 

-141.35 

482 


WKBKM 

%Error 

20.91 

-99.62 

-2.431 

76200 




KS2H 

5 



■Dl 

-144.92 

758 

-162.83 

73 

-180.46 


Exact 

-141.70 



■DB 




WM 

%Error 

20.75 







BsEgiSl 

6 

Predicted 

-169.56 

1 


743 

-162.72 

73 



Exact 

-141.38 

391 

-148.49 

1 

-140.21 

647 

-145.35 

nm 

%Error 

19.93 

-99.74 

-2.317 





HZ2H31 

-50 

3 

Predicted 

-127.98 



1 

-133.49 

1844 


■ 

Exact 

-138.53 

199 

-148.49 

1 

-139.93 

383 


73 

%Error 

- 7.62 

328.64 

8.775 

0 

-4.602 

381.46 

-3.039 


4 

Predicted 

-127.93 


-161.77 

1 




73 

Exact 

-137.79 

262 

-148.49 

1 

-138.70 

482 


73 

%Error 

- 7.16 

225.95 

8.943 

0 

-3.756 

303.112 

-3.026 

BEHl 

5 

Predicted 

-127.92 



1 

-133.49 

2028 


73 

Exact 

-137.17 



1 


567 

-145.40 

73 

%Error 

-6.74 

161.96 

8.815 

0 

-3.050 

257.67 


■Hi 

6 

Predicted 

-127.94 

855 

-161.00 

1 




73 

Exact 

-136.64 

390 

-148.49 

1 

-136.81 

647 

-145.38 

___ 73 

%Error 

-6.36 

119.23 

8.425 

0 

-2.427 

225.81 

-2.992 

miMi 
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Table A.3 Simulation results for the test patterns i.e. patterns not seen during training (adequate 

training case ) 


Temp 

Change 

(degC) 

Tran 

Time 

(min) 


Inner radius 

Outer radius 

Max Stress 

Min stress 

Max Stress 

Min stress 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 

Value 

(MPa) 

Time 

(s) 



Predicted 

-148.45 

1 

-153.16 

221 

-145.56 

73 

-148.07 

424 


3 

Exact 

-148.49 

1 

-154.46 

199 

-145.22 

73 

-148.54 

382 



% Error 

-0.02 

0 

-0.84 

11.05 

0.23 

0 

-0.31 

10.99 



Predicted 

-148.62 

1 

-153.62 

292 

-145.73 

73 

-148.70 

517 


4 

Exact 

-148.49 

1 

-154.90 

262 

-145.25 

73 

-149.28 

481 

+30 


%Error 

0.08 

0 

-0.82 

11.45 

0.33 

0 

-0.38 

7.48 


Predicted 

-148.97 

1 

-153.56 

348 

-145.61 

73 

-148.76 

582 


5 

Exact 

-148.49 

1 

-155.27 

326 

-145.26 

73 

-149.89 

566 



%Error 

0.32 

0 

-1.10 

6.74 

0.24 

0 

-0.75 

2.82 



Predicted 

-148.87 

1 

-153.74 

408 

-145.31 

73 

-149.02 

670 


6 

Exact 

-148.49 

1 

-155.90 

391 

-145.27 

73 

-150.41 

646 



%Error 

0.25 

0 

-1.38 

4.34 

0.02 

0 

-0.92 

3.71 



Predicted 


1 

-159.37 

173 

-144.96 

73 

-151.30 

352 


3 

Exact 

-148.49 

1 


199 

-145.17 

73 

-150.69 

383 



%Error 

-0.05 

0 

huh 



0 

0.40 

-8.09 



Predicted 

-148.27 

1 



-145.14 

73 

-152.97 

480 


4 

Exact 

-148.49 

1 

-159.18 

262 

-145.20 

73 

-151.92 

482 



%Error 



0.93 

-1.14 

-0.04 

0 

0.69 

-0.41 


Predicted 


1 

-161.15 

■ctoH 



-153.68 

532 


5 

Exact 

-148.49 

1 


MEM 

-145.22 

73 

-152.93 

567 



%Error 

-0.10 

0 


-7.05 



0.49 

-6.17 



Predicted 

-148.33 

1 


361 

-145.04 

73 

-154.66 

619 



Exact 


1 

-160.33 

mm 


73 

-153.81 

647 



%Error 

-0.10 

0 

0.91 

-7.67 

-0.13 

0 

0.55 

-4.32 


■ 

Predicted 

BEHEI 

1 

-160.26 



73 

-151.11 

401 


J 

Exact 

-148.49 

1 

-162.42 

199 

-145.11 

73 

-152.84 

383 



%Error 

-0.08 

0 

-1.32 



0 

-1.13 

-4.69 



Predicted 


1 

-161.75 



73 

-153.55 

503 


■ 

Exact 


1 

-163.46 

262 


73 

-154.56 

482 



%Error 

-0.14 

0 

-1.06 

4.58 

0.10 

0 

-0.65 

4.35 


Predicted 

-148.83 

1 



-145.84 

73 

-154.98 

wm 


5 

Exact 

-148.49 

1 

-164.32 

326 

-145.16 

11 

-155.98 




%Error 

0.22 

0 



0.46 

0 

-0.64 

-0.53 



Predicted 


1 

-162.72 

383 

-145.46 

73 

-155.93 



6 

Exact 

-148.49 

1 

-165.07 

391 ' 

-145.21 

73 

-157.22 

mum 



%Error 

0.20 

0 




0 

-0.82 

-2.62 

-30 






1 



wmm 

73 


3 



199 

■Eliga 

1 



-145.40 

73 





14.57 

MSEM 



11.48 

-1.05 
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WS&ffl 

Predicted 

-140.41 

255 

-148.36 

i 


479 


73 


B| 

Exact 

-142.07 



i 

-141.35 

482 

-145.38 

73 




-1.16 


-0.08 

0 

-0.67 

-0.62 

-0.05 

0 



Predicted 

-140.08 

317 

-148.35 

i 

-139.78 

549 

-145.49 

73 


5 

Exact 

-141.70 

326 

-148.49 

i 

-140.74 

567 

-145.36 

73 



%Error 

-1.98 

-2.76 

-0.09 

0 

-0.68 

-3.17 

0.08 

0 



Predicted 

-139.37 

381 

-148.38 

i 

-139.07 

628 


73 


6 

Exact 

-141.38 

391 

-148.49 

r i 


647 

5PS 

73 



%Error 

-1.80 

-2.55 

-0.07 

0 

-0.81 

-2.93 


0 



Predicted 


178 

-148.17 

1 

-141.55 

353 

-141.02 

73 


3 

Exact 

-138.53 

199 

-148.49 

1 

-139.93 

383 

-145.45 

73 



%Error 

1.80 

-10.6 

-0.21 

0 


-7.83 

-3.04 




Predicted 

-139.81 

248 

-148.30 

1 



-142.19 

73 


4 

Exact 

-137.79 

262 

-148.49 

1 

-138.70 



73 



%Error 

1.46 

-5.34 

-0.12 

0 

0.95 

-4.77 

-2.22 

0 


Predicted 

-139.01 

316 

-148.35 

1 

-139.08 

545 

-144.80 

73 


5 

Exact 

-137.17 

326 

-148.49 

1 

-137.69 

567 

-145.40 

73 



%Error 

1.34 

-3.06 

-0.09 

0 


-3.88 






-138.23 

379 


1 


625 

-145.41 

73 


6 

Exact 

EB^SI 


beeesei 

1 

EfZ&jl 

647 

-145.38 

73 


%Error 

1.16 

-2.82 

-0.04 

0 
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